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SUMMARY 


This report addresses the problem of controlling the vibrations of large space 
structures by the use of actively augmented damping devices distributed through- 
out the structure. 

When these devices introduce only a moderate amount of damping (10% - 20%), 
simple analytical formulas are obtained for either predicting the damping or de- 
signing the dampers from the application of the theory of Low Authority Controllers 
to structural control. Such control systems are shown to be always stable and not 
very sensitive to parameter changes or uncertainties when actuator and sensor 
are colocated, i,e. , each damping device is a self-contained unit attached to the 
structure and acting as an energy sink. 

A newly conceived damping device, the gyrodamper, is described and analyzed. 

It consists of a set of single gimbal control moment gyros which are actively con- 
trolled to extract the structural vibratory energy through the local rotational 
deformations of the structure. Various linear and nonlinear dynamic simulations 
of gyrodamped beams are shown, including results on self-induced vibrations due 
to sensor noise and rotor imbalance. The complete nonlinear dynamic equations 
are shown in Appendix A. 

The problem of designing and sizing a system of gyrodampers for a given structure, 
or extrapolating results for one gyrodamped structure to another is solved in terms 
of scaling laws. Novel scaling laws for gyro systems are derived, based upon 
fundamental physical principles, and various examples are given. 
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INTRODUCTION AND OVERVIEW 


In the development of large structures constructed or assembled in space , it is 
of interest to determine the effectivenss of active devices which augment 
structural damping. Natural damping in large structures can be significantly 
enhanced by incorporating into the structure energy dissipation devices which 
transform the vibratory motions into heat, either directly by mechanical friction 
(e.g. , passive dashpots) or through electrical networks (e.g. , active electro- 
dynamic actuators). In the latter case, a sensing device is required to measure 
the structural vibration velocities which are fed back to the actuator with the 
proper gain so that the colocated sensor /actuator will mimic a purely passive 
device. Such active damping devices (sometimes called augmented passive 
dampers) differ essentially from purely passive devices because they are 
effective over a broad range of structural vibration frequencies, whereas a 
passive device can only be tuned effectively for a single frequency. 

In general, the sensitivity of energy absorption characteristics of augmented 
passive dampers is far superior to what can be achieved by purely passive means 
such as mechanical friction. The advantages of augmented dampers include also 
their adaptability and reliability characteristics , primarily because electrical 
devices fare better in the harsh space environment than do viscous fluids or 
materials whose granular friction characteristics must remain constant, or at 
least predictable. 

Structural dampers can be classified to reflect the duality between translation/ 
force and rotation /torque , as well as the criterion which distinguishes where the 
reaction forces or torques are transmitted. If these reactions are transmitted 
back into the structure, the dampers are intrastructural , while if they are trans- 
mitted to inertial space, the dampers are inertM . The corresponding classifi- 
cation matrix, with some generic examples, is shown below, with additional detail 
shown in Fig. 4 of the report. 





Classification Matrix of Damping Devices and Generic Examples 

In the present study, attention is focused on a class of inertial torque-producing 
devices called gyrodampers . In their simplest form , the so-called V-CMG damper 
illustrated in Fig. 4 of this report, a pair of identical control moment gyros 
(CMGs) are mounted "back-to-back" on parallel gimbal axes so that, in this 
nominal position, their angular momenta are equal and opposite and hence cancel. 
Gimbal motion commands are also equal and opposite, resulting in symmetric 
V-scissoring motions of the CMG angular momentum vectors about a fixed axis. 
Thus, the resulting torque (vector sum of the individual rates-of-change of 
angular momentum) always lies along this axis. This torque, applied to the 
structure by the gyros, is then electronically controlled so as to mimic the effect 
of an (inertial) rotational dashpot opposing the rate-of-change of local structural 
angular rotation. • 
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A purely passive version of the V-CMG dampers were developed at Lockheed 
in 1963 for gravity- gradient attitude stabilization of the Agena vehicle and were 
subsequently flight- tested several times - very successfully. At that time, the 
V-CMG dampers were constructed from oversized rate-sensing gyros used in a 
reversed role, i.e. , that of actuators producing inertial rigid-body damping 
torques rather than sensors measuring rotational rates; see Refs. 11 and 12. 


The genesis of the gyrodamper can be understood most easily by examining 
intuitively the basic properties of the simplest inertial damper, the so-called 
proof-mass damper depicted above in the upper right-hand corner of the damper 
classification matrix. If we consider a generic flexible space structure as an 
oscillating spatial sinusoid, it is clear that the "optimum” location for a proof- 
mass damper is the anti-nodal point where translational velocities are maximum; 
these velocities are measured by the damper’s (colocated) sensor and fed back 
to an electronic control system which drives the proof-mass so as to mimic the 
effect of an inertial dashpot. As this damper is relocated closer to the nodal 
point, its effectiveness diminishes in that a higher feedback gain (i.e., higher 
damping force) is required where there is less structural motion. In the limit, 
the proof-mass damper is totally ineffective when located at the nodal point where 
there is no structural motion. 



Proof-Mass Dampers, Flywheel and V-CMG Gyrodampers 


When two proof-mass dampers are considered, symmetrically situated about a 
nodal point and acting in opposition (see diagram) they produce a damping force 
couple about the nodal point. In the limit, as they approach the nodal point, 
they can be replaced conceptually by an oscillating flywheel pivoted at the nodal 
point. The inertial "flywheel damper" represents thus the limiting case of a 
pair of very inefficiently located opposing proof-mass dampers. Indeed, not 
only are high torque levels required, but also high wheel inertia, with a corre- 
sponding high weight penalty. The flywheel damper does not appear thus as a 
very efficient device to produce inertial reaction torques about nodal points. 
However, it has an advantage at very low frequencies in that it is only velocity- 
limited whereas the proof-mass is displacement-limited. 

By contrast, the gyrodamper (located also at nodal points) can be considered 
as a "black box" torque battery, capable of producing high inertial torques 
without requiring the counter-rotation of large inertias. The essential char- 
acteristic of gyrodamper s is their angular momentum storage capability . Because 
it is not possible to store linear momentum , there is no translational analogue to 
the gyrodamper, which must therefore be regarded as a unique device for pro- 
ducing inertial torques for structural damping. 

A general property of intertial dampers is their inability to transmit structural 
vibration energy from one part of the structure to another, as can occur when 
several intrastructural member dampers are used simultaneously (Ref. 5). By 
contrast, the simultaneous use of several inertial dampers is "decoupled" in the 
sense that each one provides an independent inertial energy sink. Nevertheless, 
non-inertial dampers (other than member dampers) may be required for very 
low frequencies or global reshaping (e.g., curvature control), rather than just 
for vibration suppression. A heuristic chart, showing the range of operation of 
various types of actuators is shown below. 

> 

In this chart , potential prototype actuators are sketched roughly in the middle 
of their bandwidth (horizontal scale), and at the bottom of their vibration sup- 
pression range (vertical scale) . The devices shown do not yet exist as vibration 
control devices for space structures, and their characteristics overlap quite 
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ACTUATOR TYPE 



Range of Operation of Various Actuator Types 

significantly. This study focuses on gyrodampers because of their central range 
of operation and their potentially remarkable performance as high-torque level 
inertial rotational dampers. The principal drawback of gyrodampers, at the 
present time, is due to the nature of their principal component devices, the 
control moment gyros (CMGs). These are not off-the-shelf components, but 
rather custom-built , "few-of-a-kind", devices whose exorbitant cost includes the 
associated R&D development work. (Appendix C illustrates this point for a 
particular manufacturer who lists six CMG models, for a total of 30 units built.) 

The first section of the present study addresses structural vibrations which can 
be controlled by distributed systems of sensors and actuators with limited damping 
authority, i.e. , the control system is allowed to modify only moderately the 
natural modes and frequencies of the structure. This is typically the case when 
the usually low natural damping of a structure is augmented by passive or active 
dampers to produce 10 - 20% damping in the principal modes. While it is relatively 
easy to model such devices and analyze the characteristics of the total system, 
an a priori prediction of their effect is not feasible because no analytical 
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expression exists for the roots of high-order polynomials and for eigenvectors 
of large matrices. Thus any modification of the original structure, as occurs 
for instance when damping devices are embedded in it, will lead to a new set 
of complex eigenvalues /eigenvectors that need to be determined all over again 
by numerical processes. This a posteriori evaluation is quite cumbersome, and 
the original design decision on locating and sizing the dampers is then essentially 
based on guess work and engdr^eering judgment. 

A novel idea presented in this study is that it is possible to apply root per- 
turbation techniques to predict analytically the behavior of the total system. 

A new formula is established for the eigenvector perturbation which helps define 
the validity range of the well-known Jacobi root perturbation formula. Analytical 
formulas are then derived relating the root shifts and mode shape changes to the 
controllers* gains, and these can be used both for performance prediction and 
controller design. These topics are the principal elements of J. N. Auburn's 
Low Authority Structurual Control Theory (Ref. 1), and are presented in the 
first section of the report. 

The remaining sections of the report address the basic physical principles, 
dynamics, operational characteristics and scaling laws of gyrodampers. These 
actuators were studied primarily in the context of planar beam vibrations, using 
a 100-m GrEp beam as a preliminary design, but the extension of the concept to 
beam-like space platforms became readily apparent and is addressed in the last 
section dealing with Scaling Laws. Appendix A, prepared by D. A. Levinson, 
outlines the derivation of the equations of planar motion of a free- free beam 
equipped with an arbitrary number of single-axis V- gyrodampers, arbitrarily 
located on the beam. These equations include the nonlinear actuator dynamics 
of the gyros, and are, conceptually, a special case of more general equations of 
three-dimensional motion of an arbitrary finite-element free- free structure 
equipped with arbitrarily many gyros. (The associated computer program GDS, 
Gyrodamped Structures, is a modular "N-Body" program being developed under 
the Lockheed Independent Research Program.) 
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Appendix B on Local Rotation Modes presents a derivation of the nth modal 
rotation coefficients, i.e. , the contributions of the nth mode of a structure to 
the local instantaneous axis of rotation at a point. It is shown, based on an 
1858 theorem due to Helmholtz, that this local rotation axis is obtained as a 
linear modal combination of the curls of the individual modal deformation vector- 
fields. Finally, Appendix C illustrates the gyrodamper scaling process for a 
340 lb, 23* free-free aluminum beam whose planar oscillations are to be damped 
by one-axis V-CMG systems mounted at its free ends. 


Use of trade names or names of manufacturers in this report does not constitute 
an official endorsement of such products or manufacturers, either expressed or 
implied, by the National Aeronautics and Space Administration. 
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LIST OF SYMBOLS* 


a 

A 

A» 

B 

B 


mn 

t 

' 1 ' 


C,,C2, 


D 

\r ^ 

E 

fa(t) 

F 

F 

a 

G 

h, h 

^1 * ^2 * ” 
H 

H 

i 




subscript relative to an actuator 

dynamic matrix of a linear system (square matrix) 

dynamic matrix modified by the control system 

damping matrix (also denotes beam, with no ambiguity in context) 

element of matrix B 

scaling constants 

column matrix containing nonzero D elements (also gyrodamper 
damping gain, with no ambiguity in context) 
matrix of damping gains (or scalar for one damper only) 
element of matrix D 

orthonormal frame locally attached to beam 
Young’s modulus 

input to a distributed force actuator a 

load applied at nodal station j by distributed actuator a 

column matrix of implied forces 

elements of F 

weighting matrix 

angular momentum (and its magnitude) carried by one gyro 
angular momentum vectors of gyro 1,2,... 
distributed sensor matrix 

total angular momentum of gyro mounted on vibrating beam B 
2 

imaginary unit ( i = - 1 ) 


♦Excluding Appendix A which is self-contained; all symbols are defined therein, and 
because it involves derivation of equations, are much too numerous to be listed 
separately. 
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identity matrix, identity dyadic 
sectional inertia of a beam 
total beam + gyro inertia 
subscript relative to a nodal station 
total inertia of a gyro about its gimbal axis 
spin inertia of gyro rotor 
spring constant (proportional gain) 
gyro rate -feedback gain 
normalized gain = (K + h)/J 
gimbal torque gyro scaling parameter 
momentum gyro scaling parameter 
gimbal inertia gyro scaling parameter 
spin inertia gyro scaling parameter 


(no ambiguity in context) 


dissipated power gyro scaling parameter 
rotor strength gyro scaling parameter 
bending stiffness structural scaling parameter 
characteristic length of idealized gyro rotor radius 
beam length 

nth left eigenvector of matrix A 
mass of a V-gyrodamper unit 

subscript relative to mth mode ( no ambiguity in context) 

subscript relative to the nth mode 

mass matrix of finite -element structural modal 

total structural mass (no ambiguity in context) 

number of V-gyrodamper units 

number of actuators 

number of modes to be controlled 

number of nonzero gains D 

• ar 

number of modes 
number of sensors 
number of nodal stations 





q 

r 

R. 

R 

s 

S 

S 


n 

E 


ni 


s 

t 

T 


T 

g’ g 


"V 

T(q) 


u 


(r, t) 


X 

y 


z 

/?, 

a 

CO 

CO 

CO 

CO, 


n 


n 

1 

2 

G 


distributed actuator matrix 

peak dissipated power By gyro gimbal torque motor 

vector of modal amplitude , N components 

m 

subscript relative to a sensor 

nth ri^t eigenvector of matrix A 

electrical resistance of one winding turn 

Laplace transform variable 

matrix relating damping gains to damping ratios 

element of matrix S;S .= 6 (see below) 

’ ni ^a.n ^r^n ' ' 

strain tensor 
time 

external "passive" torque, magnitude, applied to gyro by gimbal 
motor 

total "active" torque applied to the gyro gimbal, = (h + K) 0 + T 

g 

"vehicle" torque applied by gyro system to the structure 

linear combination of system "modal torques" 

column matrix of structural deformations ( N components ) 

s 

physical vector of deformation at nodal station i ( 3 components ) 

physical time-varying vector field of deformation at point 7”, time t 

abscissa of a point on the beam 

column matrix of distributed sensor outputs 

components of y 

complex column vector 

dimensionless nth root of frequency equation for a free-free beam 

gyro gimbal angle 

vibration cyclic frequency ( rad/s ) 

nth mode cyclic frequency ( rad/s ) 

upper limit of gyrodamper bandwidth 

= k/J 

lower limit (= of gyrodamper bandwidth 
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amplitude of nth mode shape at nodal station j 

matrix of mode shapes 

matrices of generalized mode shapes 

nth complex mode shape 

damping ratio 

damping ratio of nth mode 

inertial angular rotation of a structural neighborhood about one 
axis 

instantaneous axis of rotation corresponding to deformation field 
TT (?, t) 

nth "modal rotation coefficient" (physical vector, 3 components) 
inertial angular velocity of the beam + gyro system 
spin rate of the gyro rotor 



LOW AUTHORITY STRUCTURAL CONTROL THEORY 


This theory addresses the analytical prediction of the behavior of structures 
controlled by distributed systems of sensors and actuators with limited authority, 
i.e. , the control system is allowed to modify only moderately the natural modes and 
frequencies of the structure. This is typically the case when the (usually low) natural 
damping of a structure is augmented by passive or active devices (so-called dampers) . 
For instance, starting with an undamped structure, and placing dampers in it in such 
a way that, say, 10% system damping is obtained in the principal modes, implies that 
the relative change of the corresponding complex root = iw^ is only 

10%. This change is small enough to justify first-order expansions, i.e. , to warrant 
the use of root perturbation methods . The following section will outline the basic 
formulas of the eigensystem perturbation theory. 

Jacobi's Formula Revisited and Further Generalizations 


The root perturbation formula was first derived by Jacobi (ref. 2) for infinitesi- 
mal perturbations which neglect the induced eigenvector perturbation. Formulas for 
perturbing the eigenvectors are developed here, leading to a more general form of 
Jacobi's formula. 

Consider first a square matrix A and one of its eigenvalues, say A^ , and the 
corresponding left and right eigenvectors L^ and defined by: 


A R 

= A 

R 

n 

n 

n 

T 

A 

= A 

T 

L 

n 

n 

n 


( 1 ) 

( 2 ) 
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Assuming now a change 6A in the coefficients of the matrix A sa that A 

becomes A* = A + 6A , the quantities and change correspondingly by 

and 6R such that one now obtains 
n 

(A + 6A) (R + 6R ) - (X + 6X ) (R + 6R ) 

' ' ' n n' ' n n' ' n n' 


which, after substituting (1), reduces to: 

A 6R + 6A (R + 6R ) = X 6R + 6X (R + 6R ) (3) 

n 'n n' nn n'n n' ^ ' 

T 

Multiplying the matrix equation (3) on the left by and using Eq. (2), one 
obtains the exact formula 



6A (R + 6R ) 
n ^ n n' 

(R + 6R ) 
n ' n m 



Jacobi's formula can now be obtained from (4) by assuming to be sufficiently 
small to be ignored and writing "6 — ■ d : " 


dA R 

dX = --- ^ 2. (Jacobi) (5) 

11 T -n 


The perturbation in the eigenvector itself is less obvious to obtain, due to the 

fact that an eigenvector is only defined within a multiplicative constant. Thus, some 

normalization must be introduced so that the eigendirection 6R^ can be uniquely 

defined. To accomplish this, another perturbation formula is obtained by multiplying 
T 

(3) on the left by R^ G , where is a normalization weighting matrix such that 


T 

n 


G R 


n n 


1 


' — ) 




1 
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Thus, for sufficiently small > denoted by dR^, we have 

G dR =0 (6) 

(Note: Because R^ is a complex vector, it could be normalized by dividing it by its 
length ||Rj^I| and denoted by R^ , where || denotes the usual Hermitian norm . If 
( , ) denotes then the corresponding Hermitian inner-product, it is generally true, in 

A A 

contradistinction to the real case, that (R , dR ) 0 . This motivates the use of 

' n n 

the weighting matrix G .) 

^ T 

After expanding (3), multiplying on the left by R^ G^ , and using (6), one obtains 

dA. s r”^ G dAR +r'^G AdR 
n n n n n n n 

This expression is then used to eliminate dX^ in (3), leading to: 


dR = (a - R r"^ G A - X l\ /r r'^ G - dA R (7) 

n \ nnn n/\nnn / n ' 

A second-order perturbation may then be obtained for dX^ by using (7) in (4) , 
rewritten with 6 -*■ d [ see Eq. (26) ] . 

State-Space Structural Equations and Modal Eigenvectors 

To apply the preceding formulas to the control of a structute, it is necessary 
to obtain stmctural equations in first-order form. 

In the lumped-mass finite-element appraoch, a structure is modeled by 
point-masses connected by various structural elements representing stiffnesses. 
In the consistent-mass finite- element model, rotational as well as translational 
degrees of freedom may also be taken into account at element nodes. Although 
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the results presented here are valid for this type of model , implementation of 

the consistent-mass approach is somewhat more involved and, therefore, only 

lumped-mass models will be considered , where each mass m. , at node i , is 

described by its displacement iT. . If N nodes are involved in the model , the 

1 s 

total system has 3N degrees of freedom. The dynamic equations of the 

s 

structure , which are linear for small displacements , may be written as : 

MU + Ku = F (8) 

where u is a column vector obtained by stacking the components of the N vectors 

s 

lu. , F is obtained in the same way from the local loads F. , and M and K are 

N X N matrices (N = 3 N ) . 
mm 'ms' 

The usual procedure is then to diagonalize (8) by solving the eigenproblem: 

where $ is a matrix of eigenvectors or mode shapes and co^/2.7r is the frequency 
associated with mode n . Because of reciprocity properties, K and M are symmetric 
matrices (M is diagnonal), and it is known (Ref. 3) that there exist a normalization 
of # such that: 

T 

^ M $ = I (I: identity matrix) 


and thus 


K # = [ - ] 

(Note that the elements of ^ have the dimension of the reciprocal square root of mass. ) 



3y defining the transformation: 


u = i>q (9) 

where q is a vector of modal amplitudes. Eq. (8) may be transformed into the 
classical modal equations 

[ '^1- J q = F 

By defining the state-vector: 



the modal equation is finally written in first-order form as: 



( 10 ) 


where 




"^2 1 
-^n 

N 

. _'S 

_I 

1 

O 


Complex eigenvalues and eigenvectors are now associated with the matrix A. 
(Note: the eigenvectors of A fin the usual sense) are called the right eigenvectors 
of A because the eigenprobfem is written as AR^ \ ^n * eigenvectors of 
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T 

A (in the usual sense) are called the left eigenvectors of A because the eigenproblem 

T T ■ T 

A L = A. L can also be written L A = X L . These two types of eigenvectors 
nnn n ^pnn 

eliminate the simultaneous use of A and A ,] 

2 

The eigenvalues of A are A = ±ico (i = -1 , and n = 1,2,.. .N_) and the 
o n n ’ m 

modal state eigenvectors are defined to be the right and left eigenvectors of A 
corresponding to > i*e.. 


III 

L" 

0. 

..0 

io) 

n 

o 

# 

o 

0 

o 

o 

1 

o 

o 

111 

[» 

0. 

..0 

1 

o 

o 

0 

0...0 

iio 

n 

t—i 

o 

o 


T 

T 


( 11 ) 


where the only nonzero elements occur in position n and n + N^. 

Fundamental Modal Damping Prediction Formula 

Damping may be introduced in the structure by measuring (components of) 
velocities* u^ at some sensor location r (r = 1 , 2 , . . . , and producing forces 
F^ at actuator locations a (a = 1 , 2 , . . . N^) such that, in general 


F = 
a 


-V D u 
Z-r ar r 


( 12 ) 


where the D 's are positive damping (rate-feedback) gains, collectively denoted by 

SLY 

the matrix D which is augmented by zeros to be dimensionally compatible with the 


*It should be noted again that Fj really denotes one component of the load at some 
node and uj one component of tile deformation at the same or at another node . For 
instance Fg n^ay denote the y-component of a force actuator at some node, while F 0 
is the z-component at the same node. 
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modal matrix $ , so that Eq. (12) can be written as F = -Dii , which is also equal 
to -D$q according to (9). Thus the structural Eq. (10) becomes; 

X = A'X (13) 


where 


A' 


= .A + dA 



2 


I 

' 0 


The new roots of the linear system (13) may then be obtained by considering 
the matrix 



T 

0 

0 

0 


(14) 


as a perturbation of the original matrix A and applying the formulas of Section 1.1. 

To the original (undamped) mode of frequency co^/2t corresponds the eigenvalue 
X =5 ici)^ , If the state-space matrix A is now perturbed by dA by the introduction 
of "moderate” damping (rate-feedback) according to (14) , the original mode will now 
have the eigenvalue X^ + dX^ , where dX^ is given by Jacobi’s formula (5). 

From (11) we have 


T 

L R 
n n 


2 ico 

n 


For brevity, define the matrix B = t^^n^ 
by: 


and the vector 




B = 

Z = dAR 
n 


(15) 


19 


Then from (11) and (14) we obtain: 


1 

-ico B, 

1 n kn 

for 

k = 



1 ® 

for 

k > 

N 

1 



m 


and thus 


T 

L Z = 
n 


-i B CO 
nn n 


Consequently, Jacobi's formula leads to 


dA s -i-B = D (f> <p 

n 2 nn 2 ^ ar ^an ^rn 

a , r 

where (b , <i> are the elements of the modal matrix 4> corresponding to the nth 
^an rn 

mode, at actuator nodal stations a and sensor nodal stations r. 'Since the <^'s and 
D's (positive rate-feedback gains) are real, it follows that the new root + dX^ 
has now a real part, and with the usual interpretation of damping ratio we obtain 
the fundamental modal damping prediction formixla: 


2 f CO 

^n n 


^ I) (j) <j) 

/. ar an^rn 

a, r 


(16) 


This is thus a general formula for predicting modal damping in a rate-feedback 
system of distributed sensors (r = 1,2,..., Np and actuators (a = 1,2,..., N^) . 
It confirms the intuitive idea that sensors and actuators should be located preferably 
where large modal deflections occur, i.e. , where mode shapes have large magnitude. 
Obviously, if is zero, mode #n is not observable and if is zero, it is not 



controllable. Moreover, formula (16) gives quantitative estimates of modal damping 
and stability. A very significant consequence of (16) is that when a = r , i.e. , when 
sensors and actuators are colocated , the structure is always stable since, for positive 
D's , the sum on the right-hand side is always positive for any n . 

Damper Design and System Robustness 

It is now possible to address the main topic of this project which is the design of 
a control system for damping the structure. Assuming that certain locations have 
been selected for actuators and sensors, the problem is to determine the gains D 

cLV 

which will result in desired damping ratios . First, the symmetry of Eq. (16) 
indicates that actuator a and sensor r can be exchanged without changing the value of 
^ . Thus, it is only necessary to define the D *s for, say a s r . Next, some 
D^j^'s may be constrained to be zero a priori; i.e., no feedback is intended from 
sensor r to actuator a . Thus a vector d containing all the nonzero D *s is 
constructed, each component d^ corresponding to a particular . 

By defining the matrix S of elements S^. = ^ajn^^rjn ’ may be 

rewritten as 

(fo;) = Sd (17) 

where (fw) is the vector of components (n - 1 » • • «N^) . If d contains 

components, tiie matrix S is then ^ • In general it will be possible to choose 

a priori the damping ratios of N modes if N, = N , by inverting (17): 

C Cl c 

d = S“^(JO)) (18) 

If N # N , other methods must be used involving pseudo-inverses, linear pro- 
gramming, or other optimization schemes, and this is a wide open field for investi- 
gation. In the case of colocated sensor /actuator damping units, is simply equal 

to the number of damping units, and thus the following result: in order to impose a 
damping ratio to Nq structural modes, Ng damping mits are required . 
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Imposing a priori damping ratios on the system modes of the damper-augmented 

system is not necessarily the best design procedure. As a matter of fact, preliminary 

computer simulations have shown that this approach will almost invariably lead to the 

use of "negative dampers" (some of the D„„’s turn out to be negative); i.e. , some of 

these devices have to put energy into the structure in order to achieve the desired 

damping ratios. Since such "negative gain" devices are inherently unstable (even 

though the truncated system is stable) , it is quite certain that the actual system with 

all its modes will be destabilized. Thus it is necessary to constrain the D *s by 

ar 

some of the aforementioned methods used when N is different from N, . 

c d 

This observation brings about the question of the sensitivity of the damping 
system to structural parameter changes and imcertainties, or conversely, the so-called 
"robustness" of this system. Assuming that the nth structural mode shape is only 
known within a relative error denoted by 6 corresponding relative error 

on may be obtained from Eq. (16) and is given by; 


6 (t <^ ) 
'^^n n^ 

'^n n 


Id 


2 


<p (() 
ar an ^rn 


S D. 


(f) <f) . 

a, r ar^an^rn 


6<l> 

T 


n 


n 


(19) 


In general, the sum in the numerator of (19) will be larger than the sum in the 
denominator, and a small error in the mode shapes may induce enough change in the 
J*s so as to drive the system unstable. However, when actuators and sensors are 
colocated, i.e., when D is zero if a ^ r , the sensitivity is minimum and (18) 

SiV 

becomes: 
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Complex Mode Shapes and Second-Order Formulas 

While in the undamped structure the mode shapes are real, i.e. , all the points 
of the structure move in synchrony, this ceases to be true in general when active 
control is introduced. Instead, the closed-loop mode shapes are characterized by the 
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fact that tiiere are phase shifts associated with each point of the structure. Thus the 
new mode shape may be decomposed in an in-phase part (real part) and a 90° out-of- 
phase part (imaginary part) 

^n *^n ^ 

Since u = is ttie value of u when q is equal to the nth eigensolution 

of the closed-loop equation (13) which is simply (R° + ^^n) ^ partition 

of containing the last components of R^ : 

r ' «nT 

R = Fr’ ! r1 
n L h i nj 

Thus 

= <l>^+ ^6R® (22) 

An approximate expression for the vector 6R^ may be calculated explicitly 
using Eq. (7) in which the matrix is chosen of the form 
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Defining; 
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it is found that: 
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and thus the complex mode shapes are given by 


(23) 


(24) 



(25) 


This shows that damping will in general resvilt in complex modes . Finally , 
the second-order perturbation for the roots can be obtained by using in Eq. (4) the 
5R*s given in Eq. (24) , leading to the formulas: 



m ?^n 


(26) 
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other expressions may be obtained which show the dependency of and 
upon the feedback gains . If one defines (for the a-th actuator , r-th sensor , and m-th 
and n-th modes) the following quantities (the so-called Bizet coefficients); 


CO 
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o <t> (m 5^ n) 
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n m 
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4 CO ^an rn ' ' 

n 


the new (complex) modes (for the j-th node) are related to the original (real) 
modes by 
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jn 
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and the corresponding perturbed mcxial frequencies are chatted by: 
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Structural Transfer Functions 

The modal equation of a structure containing damping devices may be written 
using Laplace transforms as: 


o 2 1 T 

s I + s* D* + pco^^ q(s) = ^ F(s) 
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and Eq. (9) is written as 


u(s) = $q(s) 

When the £’s are small, so that the conditions of the low authority controllers 
apply, the following approximating expressions can be found for the transfer functions 
relating the effect of a force component to the deflection u^ : 


u^(s)/F^(s) > lq„(s)/F^(s)l 


where 


B (p 
nm ^am 


q^(s)/F^(s) “ 
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These expressions are useful for evaluating the vibration of the structure induced 

bv onboard excitation sources. When the excitation frequency is close to one of the 

^ 2 2 
mode frequencies , then s = -co^ , and , since B^^ = 2 : 
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By comparison, the static deflection (s = 

0) is given by 
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Generalized Modes and Force Distributions 


So far the theory has been established for sensors measuring displacement (or 
velocities) and actuators applying local forces. The results may be extended to cases 
where sensors measure quantities such as local rotations , bending moments , etc . , 
and where actuators produce a distribution of loads (such as doublets provided by 
torquers , moment actuators , etc .) . 

The sensor output y^, may be defined by a linear transformation of the type 


y 


r 


2.H .u. 


J rj j 


(31) 


(For instance , in a one-dimensional structure ,* the local rotation may be expressed as: 


u , , - u - 
_ r+1 r-1 

^r Ax 

where Ax is the distance between stations r+1 and r-1.) Since u = €>q , the 
vector of sensor outputs may be written as: 

y = # q 


where 


is a matrix of generalized mode shapes. 


*In a three-dimensional structure, ’’locaj^ota^ion" mode shapes are obtained by 
"stacking" the three-dimensional curls V x ^ of the deformation mode shapes . 
See Appendix B. 


2 / 


In the same way , an actuator controlled by a single input f (t) may produce a 

d. 

distribution of loads: 


f . 
aj 


P .f (t) 
aj a'' 


(32) 


where j = 1 , . . . is a station index. (For example, a torquer may be defined by 
P . = 0 except for P , = -Ax , P . = Ax , where Ax is the distance between 

SiSi^ J. SLSi‘1 X 

station a-1 and a+1 . In this case, f (t) is the torque T (t) produced by the 

SL 3, 

torquer .) 

The generalized rate feedback control law corresponding to Eq. (12) is now 
written: 


f (t) = -SD y 
a'' ar'^r 


(33) 


and thus from (31) and (32) 

f . = -P .2 , D H . {l 

aj aj r , k ar rk K 

A R 

Defining the generalized mode shape matrices ^ and ^ by: 

p$ 

' A'*' Ti 

B = 



and the matrix B by: 
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all the formulas derived previously remain formally valid by replacing; 


by 


= 2. 
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In particular, the fundamental damping formula (16) becomes; 
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GYRODAMPERS 


The preceding chapter addressed the general theory of low authority structural 
control dealing with the damped vibrations of flexible structures controlled by distrib- 
uted systems of sensors and actuators with limited "damping authority." 

The purpose of this chapter is to examine in detail the implementation of a 
specific type of inertial rotational damper, the so-called gyrodampers , with colocated 
angular rate sensors. This is an important class of devices, whose application to the 
control of large space structures appears very promising. Gyrodampers may be 
considered as momentum exchange devices between "local structural angular momen- 
tum" and gyro momentum. This exchange process involves rate-of-change of momen- 
tum, i.e. , torque, so that these devices may also be regarded as torquers acting 
along the instantaneous axis of rotation of the deformable part of the structure to 
which they are attached. The torque applied by the gyros is electronically controlled 
so as to mimic the effect of a dashpot, i.e. , opposing the rate-of-change of local 
angular rotation. Thus the structural vibration energy is dissipated in the form of 
heat radiated by the gyro's gimbal torque motors. 

As will be shown in this chapter, the performance of these broad bandwidth 
devices provides substantial corroboration of the feasibility of achieving 10% to 20% 
structural damping with only a general approximate knowledge of the structure's modal 
characteristics. 

Basic Principles 

Consider a small rigid portion of an oscillating beam B , with an orthonormal 
frame {ej^ , > e^} attached to it as shown in Figure 1. For simplicity of illustrating 

the basic principles involved, we assume that B can only undergo planar oscillatory 




rotations about and that maintains an in ertially fixed orientation. Then the 
inertial angular velocity of B is 

% = (36) 


where d is the angle of rotation of B about e^^ . A gyro, with constant speed rotor 
having angular momentum h , is gimbaled on , the gimbal angle a being zero 
when h is parallel to -02 . The total angular momentum H of the gyro, expressed 
in the moving frame {e^^ , e^ , e^} is then given by 


H = h (sin cr e^^ - cos <x Og) + 


(37) 


where h = |hj = const. , and J is the gyro's inertia about the gimbal axis e^ . 

Let be the torque applied externally to the gyro. The equation of motion of 

O 

the gyro is then: 



(38) 


where d()/dt denotes time differentiation w.r.t. inertial space. If (*) denotes 
time differentiation in the moving frame {ej^ , 62 , e^} , and if we assume small motion 
linearized dynamics, then from (36), (37), and (38) we obtain: 


dt 


j^h(<rej^ - 62 ) + Joe^ j = h(o- 6 j^ + o'Sig x e^^ - fig x 62 )^ 


+ + Jfffig X 6g 


(39) 


“ h(aej^ - 0eg) + = T 


g 


n 


32 




r- 




r 

r 


r 


Assuming now that 
equation (39) along results in 


is a gimbal torque, projection of the last vector 


Ja = hj + Tg 


(40) 


If, furthermore, T = -ka - dn , where k > 0 , d > 0 are spring and dashpot 
constants, Eq. (40) represents a passive dissipative system driven by the input 0(t). 
This passive gyro system will therefore absorb rotational energy by opposing the 0 
motion. However, the transmission of this energy is determined by the magnitude of 
h/J , which, in most practical cases of oscillating beams or structures, is too small 
for effective damping. 

Active Augmentation of Passive Gyrodampers 

The inertial absorption of the 0 -rotational energy by the passive gyrodamper is 
limited in practice by the magnitude of h/J . To increase this energy dissipation, an 
angular rate sensor is added to the system to estimate the 0* velocity which is fed- 
back to the gyro’s gimbal motor with the proper gain so that the colocated sensor/ 
actuator pair will mimic a purely passive device. As discussed in the Low Authority 
Structural Control Theory,* this passive nature alone guarantees structural stability, 
i.e., this (idealized) device cannot excite the structure, but the structure can excite 
it. Such an actively augmented passive damper can be caricatured as a passive device 
with a "hearing aid" to magnify the 0 signal, and in this sense, is quite distinct from 
a purely active device which does not necessarily require either a colocated sensor or 
a rate feedback. 

The implementation of this concept will be illustrated for a system of two identi- 
cal, coaxially gimbaled gyros (so-called V-gyros) as shown in Fig. 2. Indeed, the 
assumption in the previous example of Fig. 1 that e^^ has an inertially fixed direction 
was only used for simplicity to idealize a simple 2 d.o.f. model using one gyro. In 
its quiescent state (i.e., 0 = 0 , ct = 0) , the one-gyro model has a stored angular 


*See Eq. (16) 
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momentum h // -62 , making the composite beam + gyro system a gyrostat, a 

generally undesirable characteristic. In addition, only infinitesimal gimbal motions 

♦ ^ 

will ensure that the body gyro torque h is aligned with the 6 -input axis e^^ for 
damped momentum exchange between beam $ -oscillations and gyro gimbal motions . 

When two identical gyros are used such that zero gimbal angles correspond to 
gyro momenta alignments h ff -e^ and h 2 / 62 (see Figure 2) , then the counter- 
rotation of the rotors cancels the internally stored momentum. In addition, because of 
symmetry, the 0 * -induced gyro gimbal motions result in the classical "V-scissor” 
pairing of the momenta h^^ and hg ; i.e., the resultant total gyro momentum (w.r.t. 
the frame {e^^ , eg , e^} ) always remains aligned with the ± e^^ beam axis. This 
allows utilization of increased gyro gimbal travel. Of course, this travel is limited 
in practice due to the nonlinear relation between gyro momentum and gimbal angle, 
i.e., 


(hj^ + h 2 > ' = 2 hsincr 


(41) 


where h = jhj^j = |h 2 | . 

For each gyro, Eq. (40) remains valid, but the gimbal torque T is no longer 

o 

realized by passive springs and dashpots, but generated electronically by gimbal 
motors and will include a feedback term proportional to the sensor output 0 . The 
gyrodynamics (see Eq. 40) of each gyrodamper are now described by: 

feedback loop 

J(f = h 0 * + (K^^ da~'^ka) = T^ (42) 


and the output ("vehicle”) torque applied to the beam by the gyros (see Eq. 41) is given 
by: 


T 


V 



= -2 h cos Off 


s -2hff 


(43) 


35 



I 


The gain k in (42) represents a weak spring ensuring a zero dc value of the 
gimbal angles cr . The classical gyroscopic torque fig x h , where j fig | = 0 » 
occurs on the gyro gimbal axes and appears as h0 in Eq. (42). (This torque would 
normally drive the gyro gimbals in a purely passive system.) 

For preliminary design purposes, the gains in the feedback term K0 - dcr in 
Eq. (42) are made large enough for the dynamics of the gyro itself to become negligible 
(i.e., the Ja inertial effect), ensuring d ~ [(K + h)/d] e , so that Eq. (43) becomes: 



D 


= -D0 

= 2h(K + h)/d 


(44) 


n 


where D can be interpreted as the damping constant of an ideal 0 -damping torquer 
attached to inertial space. This interpretation reflects the fundamental property of 
gyrodampers used as inertial energy sinks, or equivalently, as inertial dampers, in 
contradistinction to intrastructural dampers (e.g., shock absorbers) connecting two 
distant points of the beam or structure in general . 

The above considerations become more quantitative by introducing the gyro- 
damper transfer function (see Figure 3) corresponding to Eq. (42): 


where 


O -/0 = 

K' s 

c^i = 
2 

= 

G 

.2 

X = 


iK'gj 

2 2 

- CO + i CO. CO 
<jr 1 


(K + h)/j 
d/J 

k/J 

-1 




(45) 


The gyrodamper bandwidth is the interval [co^ , co^^] , where 


"o ^ 


(46) 
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Since k is a weak spring, oo^ is small. On the other hand, is chosen 
large enou^ to include the highest structural frequency which must be damped . 
Within this bandwidth. 


I 2 2, 

|0)q - O) I « COj^CO 

so that the transfer function (45) reduces to 

cr/0 s K'/w^ 

Thus Eq. (44) can now be expressed as: 


where D 

Beam Damping Simulations 

Physical and Mathematical Model. To illustrate specifically the (theoretically 
predicted) performance of an inertial rotational damper, the V-CMG damping system 
shown in Figure 4 was rigorously modeled on Lockheed's multibody hybrid coordinate 
software (ref. 4) . The model included both the noise produced by the 6 rate-sensing 
xinit (e.g., rate sensing gyros, not shown) and the effects of gyro vibration due to rotor 
mass imbalance. 

The physical model used (Figure 4) is a 10-m tubular aluminum beam, 10-cm 
diameter, 2-cm thick, whose mass is 136 kg. A horizontal V-CMG damper unit is 
clamped at its center, with the gyros mounted back-to-back in the nominal (zero 

momentum) configuration. For each gyro, the rotor inertia about the shaft is J 

2 ^ 

~ 0 . 02 1 1^-m , while the transverse inertia J (including gimbal inertia) is slightly 

over half that amoimt. The mass of each rotor is 3 kg, its diameter is ~ 15 cm, and 

its spin-rate is 400 rad/s (~ 4000 rpm) so that the angular momentum of each gyro is 

h = 8.4 N-m-s (<^6.2 ft-lb-s). Additional physical parameters are shown in Figure 4. 


-2ha = -D0 
2hK'/co^ 


(44a) 
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Figure 4. V-CMG inertial damper producing single -axis damping torque 




The mathematical multibody hybrid coordinate model used consists- of five bodies: 
The ’’damping unit box”, two identical single-gimbal gyros , and two 5-m cantilevered 
tubular beams abutted end-to-end under the damping unit. Because the motion equa- 
tions of this system constitute a precise evaluation model which must include the 
actuator dynamics, the first six (planar) cantilever beam modes were included for 
each beam, in addition to one rigid body mode (damper box rotation about the 0-axis), 
and 2 degrees of freedom for each gyro (gimbal angle a and rotor phase angle, an 
ignorable coordinate) . This 17 degree-of-freedom system was modeled open-loop in 
the nonlinear N-BODY program, with the gyro gimbal control torque and 0-input 
torque treated as external inputs. Using the built-in numerical linearization procedure, 
the constant coefficient matrices of the corresponding linearized (first-order form) 
model were then input to state-space linear systems software (LISSA 6), where the 
control loop was closed. The control equation is linear and relates the gimbal torque 
to the angular rate 6 of the device and also to the gimbal angle a and rate cr . This 
last variable is necessary to make the gyro system behave like a passive rotational 
damper . 

Briefly, in the subset of hybrid linear equations corresponding to the variables 
0 and (T, the rate-sensor feedback loop was closed as follows: 

Total System: l0=-2ha + T(q) J 

Each Gyro: Ja = h0 -i- (K0* - do-k^S (47) 

feedback loop I 

where I is total (beam + gyros) 0 -inertia, and T (q) is a linear combination of 
"system modal torques” produced by the antisymmetric vibrations of the beam. The 
q’s are the modal amplitudes associated with the beam mode shapes included in the 
simulation. (See Appendix A for a complete description of the nonlinear actuator 
dynamics for a free-free beam equipped with an arbitrary number of V-gyrodamper 
units .) 


Design and Tuning Procedure 

The first step in the design procedure is to determine the value of the 
damping gain D which will give desired performance for an ideal torquer. This 
can be obtained either from the fundamental prediction formula (16) or from 
parametric simulation studies of the actual system. Using the latter, a damping 

4 

gain D = 1*2 10 N*m's was found desirable. To implement this with the gyro- 
damper described earlier (h = 8.4 N*m, J = 0.01 kg m^) with a desired band- 
width = 1400 rad /sec (about twice the frequency of the highest mode to be 
controlled), Eqs. (44) and (45) are used to determine the two control gains: 

d = 

K = 

Without this feedback loop, the equivalent damping gain of the passive system 
would have been D = 10.1 N-m-s. 

Numerical Simulation. The automated modeling procedure is a generalization of 
the ”two-stage eigenanalysis" method used by Draper Laboratory (ref. 5) in that the 
original model includes not only structural modes, but also the dynamics of other 
components of the system (e.g., interconnected bodies, gyros, etc.). Of the 34 
(complex) eigenmodes corresponding to the original 17 degrees of freedom physical 
system, 12 system modes are essentially the first 12 classical free/free-beam modes. 
Because the damper imit is located at the exact center of the 10-m beam, the six 
symmetric modes are undamped. Of the remaining six antisymmetric modes, only 
the first three are physically meaningful because of bandwidth limitations (~ 100 Hz) 
inherent to rate-sensing gyros needed for the K0 feedback term in Eq. (47). The 
corresponding system modal frequencies, damping coefficients and decay time are 
listed in Table 1. 


= 14 N'm-s 

Dd i, ~ 1 XT 

- h = 10 N.m.s 


(48) 
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TABLE 1. SYSTEM MODAL DAMPING 


Frequency (Hz) 

Damping (%) 

Decay Time (s) 

19.9 

15.09 

0.053 

57.8 

5.35 

0.051 

113.9 

2.93 

0.048 


This system was simulated for 1 s , the initial conditions being a 1-mm deflection of 
the beam tips in opposite directions. This static deflection was decomposed in the 
original cantilever modal basis to obtain the initial conditions (0) for the beam 
modal amplitudes in the closed-loop hybrid linear system . The ensuing damped 
oscillations of the tip-deflection and its modal components <lj^(t) are shown in 
Figure 5 for the first 0.5 s of simulation. Without sensor or actuator noise, tip 
deflections diminish to a few millimicrons at 1 s, an interesting but not physically 
meaningful result. In practice, sensor and actuator noise drive the gyro damper as 
an active excitation source. However, the resulting vibrations turn out to be below 
the level required by precision space structures. (For optical or rf systems, 
this level varies between 1/10 to 1/50 of the wavelength.) 

With rate-sensor noise (obtained from bandpass filtered white noise) having 
mean amplitude 0.5 /urad/s (an order of magnitude higher than for the best rate- 
sensing gyros) , the residual tip-deflections between 0.9s and 1 s average a few 
hundredths of a micron. This noise, and its effect on tip-deflection, gyro gimbal 
torque, and spurious 0 -damping torque are shown in Figure 6. For actuator noise 

due to gyro rotor vibration, one compliant bearing was assumed for each gyro to 

“8 

allow the coning motions of the rotor shaft due to a typical l-/nn. (~ 2.5 10 m) mass- 

center offset of the rotor. The resulting level of spurious damping torque 
-4 

(~ 16 10 N-m) is about half the peak amplitude of the spurious component due to 
sensor noise (Figure 6) . Finally, the time-histories 0(t) and a(t) are shown in 
Figure 7 . These results are scalable (using larger gyros) to different vibration levels . 
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Figure 7. Local beam rotation and gyro gimbal angle 
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Gyrodamping of a 100-m GrEp Beam 


The simulation studies described in the previous subsection established feasi- 
bility of the V-gyrodamping of a relatively massive, thick and stiff aluminum 10-m 
tube, including rate sensor noise and gyro-rotor imbalance effects . The damping 
unit was placed at the exact center of the beam, thus leaving (by definition) the 
symmetric modes undamped. In addition, Lockheed's general-purpose multibody 
hybrid coordinate software (in the sense of Likins) was used, requiring cantilever 
mode shapes for each half of the beam "clamped" under the centrally located 
V-gyrodamper unit. 

As can be inferred from the previous example, the numerical methodology 
employed is relatively complex. This follows from the use of generalized multipurpose 
software constrained to use cantilever modes for flexible structures attached (or 
gimbaled) as terminal "appendages" to arbitrary topological trees of interconnected 
rigid bodies. With this software, a flexible beam with a (relatively) small -mass 
V-gyrodamper attached at one end can still be modeled, but an excessive number of 
higher order cantilever modes are required to converge to a physically useful set of 
free-free modes, as was done in the first example of the massive tubular aluminum 
beam. This numerical phenomenology can be loosely identified as the pathology of 
the "tail wagging the dog." Because of this phenomenology (in addition to numerous 
other reasons derived from the rapidly evolving design changes associated with current 
and future large space structures) , development of a new generation of "multibody" 
software was initiated xmder the Lockheed Independent Research program these last 
few years . The basic idea involved is the "turning inside-out" of the original hybrid 
coordinate topological tree model originally developed in 1965 by Hooker and Margulies, 
and subsequently extensively studied and documented by Likins and many others 
(Ref. 6). 

Briefly stated, the "inside-out" model consists, in the first phase, of a single 
free-free flexible structure, having large-angle rigid-body modes, and containing an 
arbitrary number of arbitrarily gimbaled rigid bodies to model reaction wheels, gyros. 
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gimbaled variable-speed wheels, etc. (In the second phase, additional gimbaled 
flexible appendages will be attached to the main flexible body .) Because of the mathe- 
matical complexity of these models, only equations of motion for a simple 
"submodel" are given in Appendix A . This submodel is limited to the planar 
vibrations of a free-free flexible beam equipped with an arbitrary number of arbitrarily 
located V-gyrodamper units . 

To illustrate this unique simulation capability, a second example of a large 
100-m graphite epoxy beam equipped with one and two V-gyrodamper units was 
considered, including comparisons between damper locations, low authority structural 
control, and "ideal torquer" simulations, and the effects of nonlinear gyro actuator 
dynamics. The physical parameters of the beam were chosen as follows: 


Length 

Outside Radius 
Wall Thickness 
Young’s Modulus 
Density 
Mass 

Sectional Inertia 


L = 100 m 
r = 10. 55 cm 
e = 2. 275 mm 
E = 3.45 10^^ N/m^ 
p = 1607 kg/m^ 

M = 239.7 kg 
I = 8. 125 10"® m^ 


Kg = (El) = 5.973 10"^ kg ^^^-m ^^^-s 

The V-gyrodampers used on the above beam were chosen to have the following 
(conservative) characteristics (for each V-gyrodamper unit): 


Total damper unit mass 

m 

= 8 kg 

Gimbal Inertia 

J 

= 0.021 kg m' 

Rotor Spin Rate 


- 400 rad/s 

Rotor Angular Momentum 

h 

= 8.4 N-m-s 

Maximum Gimbal Torque 

T 

G 

= 3.2 N-m 

Peak Power 

P 

n 

= 8W 

Maximum Gimbal Angle 

O’ 

max 

= 1 rad 


Maximum Gimbal Rate 
Maximum Gimbal Acceleration 

Feedback Gains ^ 


a = 2.5 rad/s 

max „ 

a =40 rad/s 

max . 

K = 10 N-m-s/rad 


d =6.27 N-m-s/rad 
In the first simulation example of the aluminum tubular beam, time histories 
were displayed for all the physically relevant system variables, i.e. , 


• Beam tip deflection 

• First three modal amplitudes 

• Rate sensor noise output 

• Noise-induced gyro gimbal torque 

• Noise-induced beam damping torque (including effects of gyro rotor 
imbalance) 

• Noise-induced beam tip deflection 

• Rigid body mode beam rotation 

• Gyro gimbal angle travel 


Because of the multivariate aspect of the problem, results for the GrEp beam 
will be limited* to time-histories of the modal amplitudes of the first two free-free 
beam modes (first symmetric q^^ , first skew- symmetric ) . From these modal 
amplitudes one can reconstruct either the beam deflections ^ local 

rotations 0 = . (See Appendix B.) For all the simulations shown, the initial 

conditions for the q*s were obtained by decomposing a 1-cm beam tip deflection into 
the free-free beam modal basis . 

Figure 8 (b), (c), (d) shows the comparison between nonlinear actuator dynamics, 
low-authority structural control algorithms, and the ideal torquer model interpreted 
as the application of a damping torque directly against inertial space. (For practical 
purposes, this model corresponds essentially to linear actuator dynamics with infinite 
bandwidth.) There is a perfect match between this model and the low-authority control 
algorithm [Figure 8 (c) and (d)] , while the "modulation" occurring around 1.75 s 


and d were obtained as in (48), for the new value of J. 

*The remaining system variables have essentially the same qualitative behavior as in 
the first aluminum beam example. 
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Figure 8, Gyrodamping of a 100-m GrEp beam — comparison between nonlinear actuator dynamics, 
low authority structural control, and ideal torquer dynamics. 
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for the nonlinear actuator dynamics corresponds to a momentary gyro momentum 
saturation. Figure 8 (a) shows the effect of- nonlinear actuator dynamics on the first 
mode, which exhibits approximately 19% damping, in accordance with the theoretically 
predicted values as shown in Table 2. The difference between predicted and 
actual values is due to the relatively high level of required damping. 


TABLE 2. SYSTEM MODAL DAMPING 


System 
Mode No. 

Freq. (Hz) 

Predicted 
Damping (%)'^^ 

Actual 

Damping 

Actual Decay 
Time (s) 

1 

0.405 

19.66 

19.2 

2.04 

2 

1.14 

19.66 

13.6 

1.03 

3 

2.05 

19.66 

8.7 

0.89 

4 

3.04 

19.66 

6.0 

1.14 


(a) Obtained from Eq. (65). 


Figure 9 shows the comparison between V-gyrodamper locations when two 
units are used, via their effect on the first two modal amplitudes. As can be intuitively 
expected, the centrally located unit has essentially no effect on the first symmetric 
mode [compare Figure 8 (a) with Figure 9 (a)] , and some effect on the first skew- 
symmetric mode [compare Figure 8 (b) with Figure 9 (b)] . The smallness of this 
effect reflects the considerably smaller strain energy stored in this mode as compared 
to the first symmetric mode for the chosen initial condition of 1-cm beam tip deflec- 
tion. Also, because the introduction of damping into the system couples the modal 
amplitudes , the transients in all the simulations shown are essentially dominated by 
the first symmetric mode. Finally, Figures 9 (c) and 9 (d) show the radical improve- 
ment when the two gyrodampers units are placed one at each end of the beam. This 
again is in agreement with physical intuition in that the maximum slope changes will 
occur at the ends of a free-free beam . 

Gyrodampers for 3-D Structures 

The V-gyrodampers discussed in the previous sections are an example of 
gyroscopic single-axis damping of planar beam vibrations for which the^e units were 



designed. If transverse beam vibrations were also present, the gyros would keep 
their individual momenta antiparallel (total zero momentum) , and therefore would not 
participate in transverse beam vibrations except through a weak gimbal inertia effect. 
For two- or three-dimensional vibrations, different gyrodamping configurations are 
required, and many combinations are possible. In particular, two single-axis 
V-gyrodampers, clamped on the same beam with their gimbal axes at right angles, 
would extend the damping procedure for planar vibrations to both planar and transverse 
beam vibrations in a simple, essentially decoupled way . 

For two-dimensional beam vibrations, the damping torques which must be 
applied by the gyros will span a plane orthogonal to the beam . The two orthogonal 
V-gyrodamper units mentioned above are thus a particular implementation of a two- 
axis torque control system. There exist many other gyro configurations for two-axis 
control. For instance, two V-gyrodampers may be arranged with their gimbal axes 
parallel (i.e. , orthogonal to the control plane) and such that the "scissoring*' motion 
of each pair occurs along perpendicular axes in the control plane, as shown in Fig. 10a. 
An important quantitative characteristic of any gyro system is its momentum envelope 
(ref. 7) , i.e. , the set of all its achievable angular momenta. For the configuration 
in Fig . 10a , this momentum envelope is a square of side 4h , provided that ± 90 ® 
gimbal travel is allowed for all gyros . However , if the four gyros are not paired but 
made independent, a larger momentum envelope is obtained, i.e. , a disk of radius 
4h, as shown in Fig. 10b. The increase (by a factor tt) of the momentum envelope 
area must be traded off against the increased complexity of the damper actuation 
control laws . Finally , it can be shown that three independent coaxial gyros (whose 
momentum envelope is a disk of radius 3h) is sufficient to achieve 2-axis control. 
However, the requirement of unlimited gimbal travel in conjunction with bandwidth 
limitations may lead to unfeasible implementations of such configurations whose main 
advantage is reduced damper mass (three gyros instead of four , grouped in two 
V-pairs). These are open questions at present, and would require further in-depth 
study. 

When the most general case of a large space structure is considered , deforma- 
tions of a local material neighborhood are such that the instantaneous rotation vector 
may take any direction in space. This is the case, for example, when planar shearing 
modes of a platform are considered. The instantaneous rotation vector is then normal 
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Figure 10. Momentum envelopes for 2 V-gyros and 4 single gyros 
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to the platform , while for bending or twisting modes it is essentially tangent to the 
platform. A gyrodamper configuration for this application must therefore be capable 
of achieving 3-axis (damping) torque control , and hence its momentum envelope must 
be a three-dimensional region. Here again, three pairs of orthogonally mounted 
V-gyrodampers may offer the simplest implementation of controller actuation, with 
a momentum envelope being a cube of side 4h. If these units can be mounted on the 
edges or corners of a large platform, within easy reach of the Shuttle for replacement, 
the failure-reliability problem is alleviated. If not, failure of any single gyro auto- 
matically disables the V-pair to which it belongs, resulting in loss of control of one 
axis . 

A general approach to using six gyros for three-axis control is to abandon the 
concept of V-pairs, and moimt them as independent gyros whose gimbal axes are 
equidistributed along the generators of a cone, as shown in Figure 11. The corre- 
sponding momentum envelopes are shown in Figure 12 for the case of zero, one, and 
two gyro failures .* With no failures , and for a 30 ° cone half-angle , the momentum 
envelope is an oblate spheroid of equatorial diameter 11.21 h and polar diameter of 
6 h. In addition to the increased momentum capability of this configuration, these 
conical systems of gyros , none of which play a preferred role , can be made inherently 
tolerant to successive device failures, with minimum loss in system performance. 
However , special control laws are required to achieve this , and a tradeoff between 
complexity of implementation and structural bandwidth constraints must be analyzed. 

In the context of attitude control , on-going studies funded by the Lockheed Independent 
Reasearch Program have documented the formulation, development, and 
design of nonlinear optimal control laws for arbitrary systems of identical 
single-gpmbal control moment gyros. There appears to be a direct applica- 
bility of this technology to NASA’s interest in the control of large space structures, 
in particular when the attitude control and structural control problems must be 
synthesized. While these problems are outside the scope of the present study, it is 
quite conceivable that dual-purpose gyro systems can be designed to solve the attitude 
and structural control synthesis problem. 

*These envelopes are boundaries of almost convex regions; inside the white holes, 
the surface folds inward slightly, resembling a shallow depression. Since the 
surface is represented by a curvilinear parameter net , and since this net becomes 
singular inside the white holes, inside details are not show. For further details 
on this point, see Ref. 7. 
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Figure 12. Two views of the angular momentum envelopes of the 6-CMG cluster 
(corresponding to conical array of 6 CMGs shown in Figure 11) 
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SCALING LAWS 

The purpose of this chapter is to define the basic parameters that determine the 
performance of a structure controlled by gyrodampers, and obtain formulas enabling 
the designer to extrapolate the results obtained with a particular structure/gyros 
combination to a structure of different characteristics . 

Scaling Parameters for Control Moment Gyros 

Control moment gyro system performances are limited by four basic constraints: 

(1) The gimbal angle a : limited because of the nonlinear characteristic of the 
gyros. For instance, for V-gyros, the output torque is proportional to 

2 h cos O’ ; thus for a near 90°, the system becomes very ineffective. 

(2) The gimbal angular rate a : limited because the output torque (sho) has to 
be transmitted from the wheel to the shaft and is thus only as good as 
whatever structural part connects the rim to the shaft. 

(3) The gimbal angular acceleration u : limited because the maximum torque 
available from the gimbal motor has to overcome the inertia about the 
gimbal axis of the gimbal /wheel/spin-motor combination. 

(4) The dissipated power Pp : limited by the heat dissipation capabilities of 
the gimbal torque motor. 

Given these four constraints, the optimum gyro will be. the one with the greatest 
possible momentum h . On the other hand, the mass of the total damper device 
should be minimum. Therefore, the following problem must be addressed: assuming 
that for a given mass m , some sort of "optimal" design has been achieved, how is 
the performance of the system affected by choosing a different mass ? The materials 
used being assumed the same, the problem is equivalent to choosing a different 
characteristic length, i , of the system. Under these conditions of pure geometric 
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scaling, it will be shown that there exist some invariant quantities, called henceforth 
the gyro scaling parameters (GSP), which relate some fractional powers of the mass 
to the gyro characteristics and performance. 

Very simple physical considerations need to be introduced, bearing in mind that, 
under this geometric scaling, the mass of any one part stays proportional to the total 

o 

mass m. First, the inertia J of the wheel about its spin axis is proportional to mjt , 
and since m varies with i 

Js " (49) 

(The symbols » ^2 * designate constants invariant under geometric 

scaling.) The same type of relation holds also for the total inertia J about the 
gimbal axis. 

Next, calling the spin rate, the centrifugal force on the rim is ~ mfi f . 

® 2 2 2 s 

Since material strength is ~ £ , it is necessary that < C„ . 

S 2 

Since h = 0 , this may be written as: 

s s 

and using (49) shows that the maximum achievable momentum is 

h = (50) 

Now the output torque, ho , being limited by the wheel/shaft/bearing strength, is 
such that 


ho < 
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Combining this with (50) , an expression for the maximum gimbal rate is- 


°^MAX = ^5”^ 


-2/3 


(51) 


I 


Some consideration is now given to the electrodynamic actuator. In a classical 
type, it is assumed here that the force produced is proportional to V/R^ , where 
V is the applied voltage, the resistance of one winding of length 2 vrd . Since 
= p 2 7 rd/s ~ l/£ , where p is the resistivity of the conductor and s its cross- 
sectional area, the torque produced on the gimbal is then 

Tq = Cgm^/^ (52) 

(If the current, instead of the voltage, was maintained constant, Tp, would then be 
1/3 

proportional to m .) Finally, the maximum power is determined by the heat 
dissipation which would be related to radiating surfaces , and thus 

Pp - (53) 


) 


i 

I 


It is convenient to express Eqs. (49) to (53) in terms of the following gyro scaling 


parameters: 





J-1 

o 

5/3 
m ' 




o 





= J-1 

5/3 
m ' 




= h 


Gyro Scaling , 


H 



Parameters (GSP) ^ 



°"max 

m2/3 




= 

-2/3 
m ' 



G 

G 





= P 

-2/3 

m 



P 

P 




SI Units 


(spin inertia) 

^ 2/3 -2 

Kg m 

■ i 

(gimbal inertia) 

^ 213 -2 

Kg m 

: 1 

(momentum) 

^ -1/3 2 -7 

Kg m sec 

(54) 

i 

(rotor strength) 

„ 2/3 -7 

Kg sec 


(gimbal torque) 

^ 1/3 2 -2 

Kg m sec 


(dissipated power) 

^ 1/3 2 -3 

Kg m sec 

1 
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other characteristics may be derived from these expressions, such as the spin 

rate = h/J , the maximum gimbal angular acceleration o,,. ^ ~ T^/J , or the 
S S JVi Cj 

output torque Ty = h : 


a = 


K„K^ m 


-1/3 


or 


MAX 


K^Kjm 


-1 


T = 

-'•TT 


K^K^m 


2/3 


(55) 


Table 3 gives the values, in SI units (see (54)), of the GSP’s of a classical 
Sperry control moment gyro*, and of the smaller V-gyrodamper^ used in the beam 
simulations of Section 2. 


TABLE 3. SPEBRY CMG AND GYRODAMPER SCALING PARAMETERS 



m 


Kj 

Kh 

Kr 

Kg 

Kp 

Sperry 

78.75 

673 

861 

2.0 

36.7 

.43 

.87 

Gyrodamper 

8 

1524 

1524 

.53 

8 

.8 

2 


The last five GSP's are directly related to the performance of the gyro control 
system: the larger they are, the better the performance. Once the design has 
maximized these parameters, then the mass m may be chosen to fit a particular 
application. 


♦Calculated from manufacturer's data (verbally transmitted). 
fCalculated from Lockheed in-house preliminary design (see Fig. 4). 
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V-Gyrodampers Scaling Laws 


The performance of the V-gyrodamper (gyro + feedback loop combination) will 
be now evaluated. The four conditions or limits described previously are first formulated 
in terms of the GSP using (54) and (55): 


''max 

2) |i| < 

3) |JI s KgKjm-' 

4) P s K 

' P P 


(56) 


Then the relationship between the amplitudes of the local rotation angle 0 (due 
to structural vibration) and the gimbal angle cr is obtained from the gyrodamper 
transfer function given in Eq. (45) : 


0(OJ) 



iK'co 

2 

O) + w 


(45) 


Within the bandwidth of the system, this transfer function is essentially constant 
and the ratio of the amplitudes of a and 0 is then; 



As was shown before, Eq. (47), this system acts as a rate -feedback controller 
since its torque output is 


= -D6 


where 


Thus 


D = 2hK'/o>, 


a/e = D/2h 


(57) 


I a I by ao) 


Replacir^ in (56) <1*^1 


T. u T 2 3 

P by Ja o) 

P 


and using (57) and (54) , the conditions (56) become: 


gimbal travel 
gimbal rate 
gimbal acceleration 
power 


0 < 2 D Kjj m / 

e < 2 


0 < 2 D 
0 < 2 D 


K„ K K m^'^^ oo~^ 

H VJT 


1/2 1/2 5/6 -3/2 

Kj m 0 ) 


(58) 


This indicates that the magnitude of the initial vibration present in the structure 
(and reflected by the local rotation angle 0) that can be handled by the gyrodanipers is 
bounded by four limits, three of which are frequency dependent as shown in Figure 13. 
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Figure 13. Typical region of operation of gyrodampers 



More amplitude can be tolerated if the damping gain D is decreased, but then, 
because of the low authority controller formula (16) , the damping ratio will also 
decrease. Conversely, if more damping is needed, D must be increased, but then 
less vibratory amplitude can be removed. The quantity ^0 appears therefore 
naturally as a better variable to express the conditions (58) . Since f is not only 
related to D but also to the properties of the structure by Eq. (16), and because the 
structure will vibrate mainly at frequencies near its modal frequencies, some scaling 
parameters will first be derived for the structure itself. 

Structure Scaling Parameters 

Considering the case of a uniform beam damped about one axis by one 
V-gyrodamper (as was the case for the previous simulations), the generalized 
formula (35) leads to: 


»2 


2 t 00 = Dd> 

m n ^n 


(59) 


where represents the nth "rotation mode shape" at the actuator (and sensor) 
location. (The indices a and r have been dropped to simplify the notation .) In 
this case 



9<?>„ (x) 

ax 


(60) 


where <p^ (x) is the usual deformation mode shape at a distance x from the origin of 
the beam. Assuming that the gyrodamper is at the origin of a free-free beam, 
classical beam theory shows that* (ref. 8): 




( 61 ) 


*In the reference, and L are respectively called and f . 
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where /3^ is the dimensionless nth root of the frequency equation and depends only 
upon the boundary conditions, L is the beam length, M its total mass, I the cross- 
sectional inertia and E the Young's modulus. Defining: 


Kg = 


Eq. (59) may be expressed as: 

J = 2D K 

s 


(62) 


(63) 


Thus this particular structure can be characterized in terms of its length and mass, 
its bending stiffness parameter and its dimensionless frequencies . 

V-Gyrodamped Structure Scaling Laws 

Before assembling the formulas of Eqs. (58)— (63) for obtaining the scaling laws 
for the total system, two new parameters have to be introduced. First, for perform- 
ance and reliability reasons, there will generally be more than one damper controlling 
a particular direction in the structure. Thus let N be the total number of V-gyro- 
dampers. Second, a reasonable mass ratio should be considered between the damper 
system and the structure itself. Let p be this ratio: 


4 Nm 


(64) 


We note that with N dampers, the damping ratio will be N times greater ; thus 
(63) becomes: 


f = 2DNK M 
^n s 


(65) 


Using (64) and (65) the following scaling laws are obtained from conditions (58): 
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( 66 ) 


1/3. -1/2 -3/2 
L co^ 


Gyro parameters Mass Structure parameters 

ratio 


Here 0^ represents the angular vibration at the modal frequency o)^ and 
the corresponding damping ratio introduced by the gyrodampers . These damping 
ratios are independent of the frequency as shown in (65) and this is a particular 
property resulting from using damping torques, which involve local rotation modes, 
and from locating the gyros at the extremities of the beam. If the gyros were situated 
differently, then Eq. (65) will be frequency dependent. For instance, at the center of 
the beam, = 0 foi" all the symmetric modes. 

The conditions (66) may be plotted in the same way as displayed in Fig. 12, 
except that now the quantity 0^ is plotted instead of 0 . Since the modal frequen- 
cies are related to the structure parameters by Eq. (61), it is possible to reformulate 
Eq. (66) as; 
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These scaling laws are valid for beams and beam-like structures. - For instance, 
they could be applied to a rectangular platform and its bending modes controlled by 
gyrodampers situated along the edges or corners. There are of course many ways 
in which to use these laws. Some particular applications (see Figures 14 and 15) 
are shown below for the following values of the parameters : 



■St 

■^G 

■'j 

■'p 

K 

s 

2.5 

23.4 

1.03 

800 

2.05 

1.5 10"® 


The above parameters correspond to the following CMG design: m = 8 kg, 

J = 0.04 kg-m^, h = 40 Nm-sec, 6- = 5.85 rad/sec, T^ = 4.12 Nm, = 8 watts. 

° max u p 

The nominal structural design is chosen as p = 0.1, M = 800 kg, L = 200 m, and 
the first two bending-mode roots are: 

/3^ = 4.73 , ^2 = 

Sometimes it is useful to use the parameter , where is the decay 

time of the damped nth mode, instead of 0^ , when mission requirements impose 

that the structure be quiet after a certain time. Since ’ foj^^^ulas 

similar to (67) may be obtained by using Eq. (61) to eliminate . Figures 16 and 
17 show the influence of the length L and of the number N of dampers or the param- 
eter ’ This last scaling is particularly interesting: it shows that, for damping 

purposes, it is more advantageous to use several little gyros than a large one. This 
is due to the fact that the performance is limited by the rotor strength (a ) and 
the gimbal motor capability (a ) . This is in contradistinction to the use of gyros 

for attitude control where the momentum capability (i.e. , o ) or the power consump- 
tion may be dominant factors, in which case the curves of Figure 17 show that one unit 
is preferable to many. However, this last condition becomes important as N increases 
so that there is, in general, an optimum number of dampers for a given structure. 






(2) NOMOTHETIC HOMOGENEOUS STRUCTURE (E = CONSTANT) 

(M ~ , Kj ~ L'^) 
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HO-MO-THET'IC, adj, homothetic fig- 
ures. Figures so related that lines joining 
corresponding points pass through a point 
and are divided in a constant ratio by this 
point. 

homothetic transformation. See simili- 
tude — transformation of similitude. 


68 


(WHERE THE C's ARE CONSTANTS INDEPENDENT OF L) 

Figure 15 . Particularized scaling laws — homothetic homogeneous structure 








CONCLUSION AND RECOMMENDATION 


The feasibility of controlling the vibrations of a structure by colocated sensor/ 
actuator damping devices has been demonstrated in this study using beam or beam -like 
structures and gyrodampers as an example. Scaling laws have been derived using the 
Low Authority Structural Control theory and the physical properties of control moment 
gyros. These laws provide excellent tools for the design of low authority control sys- 
tems for beams or beam -like structures and permit extrapolation of the results for 
structures of different characteristics. 

The control of structural flexiblility implies, in some sense, the ability to 
construct, in space, large space strucutres whose behavior is determined pre- 
dominantly not by passive characteristics but by active control systems. In this 
preliminary investigation, it was useful to direct consideration to a few simple 
models in the hope that generic characteristics of the control process will emerge 
and be directly applicable to more detailed, realistic systems. Experience with 
detailed structural models, however, has already indicated that system per- 
formance metrics, such as rms surface shape, rf or optical performance mea- 
sures, modal damping, etc., can be extremely sensitive to changes in the 
structural model, actuator and sensor placement, and sensor measurement error. 
These sensitivities may also result from the detailed way a portion or all of the 
structure is built, the limitations of sensor measurement, or the disturbances 
associated with specific actuator implementations. Therefore, while consideration 
of quite simple structures is a necessary first step in obtaining knowledge of 
fundamentals, it should not generally be expected that control strategies 
developed in this way can be immediately implemented on real, high-detail 
structures. For the latter, questions of local versus global control, sensor/ 
actuator placement strategy and error, model accuracy and identification, modal 
density, and computational requirements will be dominant considerations. 
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The positive results obtained in this study strongly motivate further investigations 
on more complicated structures, since the fundamental phenomena are now better 
understood. In particular, a rectangular reticulated space platform made of the 
102-in. NASA graphite epoxy columns will provide a more generic example on which 
to evaluate both transient effects and external disturbances, while being controlled by 
gyrodampers. 

Also an investigation into other kinds of damping devices is needed, in particular 
into the linear-inertial type ("proof mass" damper), so that tradeoff studies between 
various actuator types could be available. 

Finally, the analytical and simulation results obtained in the present study will 
serve as a basis for design of a laboratory experiment which will attempt to demon- 
strate the practical implementation of structural control . 
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Appendix A 

EQUATIONS OF MOTION FOR A BEAM CARRYING V-GYRO UNITS* 


The differential equations governing the motion in an inertially fixed plane A of a 
beam B equipped with an arbitrary number of V-gyro units can be derived as follows: 
First, to put B into a general configuration in A, introduce A^^ and A^ (see Fig. A-la) 
as perpendicular lines fixed in A, intersecting at point O . Align the neutral axis v 
of B with Aj^ such that the intersection of v with one end of B lies at O, Next, 
treating B as if it were rigid, move this point of B to a point Q of A (see Fig. A-lb) 
The position vector p^'^^ of Q relative to O is then given by 

pQ/O ^ X a^ + Y ag (A. 1) 

where a^^ and vectors parallel to A^^ and A 2 as shown. Now rotate B 

an amount 9 relative to A^ , letting and be lines respectively parallel and 
perpendicular to v and intersecting at Q . Finally, subject B t(^a flexural deforma- 
tion relative to (see Fig. A-lc), so that the position vector p from Q to a 
point P of can be expressed as 

P /o 

P = X Cj^ + y c^ (A. 2) 

where c. and c are unit vectors respectively parallel to C. and C , as shown. 

A is presumed to be a principal plane of flexure for B. 


*ByDavidA. Levinson. 
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Next, it is assumed that y can be written 




i=l 


where . . . , are functions of only the time t , and (p^ , 
upon X. Defining the functions as 


(A. 3) 


, . . , <f}^ depend only 


Uj^ ^ X cos 0 + Y sin 0 


U2 A -X sin 0 + Y cos 0 


(A. 4) 
(A, 5) 


U3 4 e 


(A. 6) 


Ug Aq^ (i n) 


(A. 7) 


A P A P ~ 

one then can express the velocity v and acceleration a of P in A as 


A P 

V 


= lu, - U3 




(A. 8) 


and 


A P 
a 


( n n 

U1-U3 2'^3-Hi h- 


U2U3-U^x Cj 


|u2+ I Ug^. ^ + U3X + U3U^-u2 I q.0.jc 


(A. 9) 
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A P A P 

Partial velocities » • • . » ys+j^ associated respectively with , 
then can be determined from Eq. (A. 8) and are found to be (Ref. 9, p. 26) 




A P 

5^1 ° 5i 


j = l 




Letting L and p stand for the length and the mass per unit length of B , one can 
form the contribution ^F| j to the generalized inertia force associated with 

(f = l, ...,3 + n) from the relation (Ref. 9, p. 89) 


("J), 


A P A P - 
• a p dx 


Equations (A. 9) through (A. 14) give 


n n 
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-UgM 


(A. 16) 


n 

(f|)^ = -u,m- J 


n 




Ugd-UgU^M+U; 


i = l 


2 “i*i 


i = l 


n n n n n 

ini ' 2 %»i - 2 2 2 2 ^ij 


i=l i=lj=l 

n n 


i=l j=l 

- ''2^3 2 “1*1 ■ "2'* ' 2 "3+1 «I - “3' - "3"l‘> 


i = l 


i = l 


n 


where 




i = l 


n 


+ U" 


y q. E.. 

Z, 1 1] 


3+i^iJ-«3Nj-''3«l*J 


(j = 1, ... , n) 


i = l 


M A J p 6x 


d A J X p dx 


I A J p 6x 


j-j 

A / pdx (i = 1, ... , n) 


(A. 17) 


(A. 18) 


(A. 19) 


(A. 20) 


(A. 21) 


(A. 22) 
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n 


n 


N 


i # Jx<^. pdx (i = l,...,n) 


(A. 23) 


XJ 

Ey A f (t>^(j)^p6x (i, j = 1 n) 


(A. 24) 


n 


If one chooses for (i = 1, n) the so-called free-free modes of a uniform 
beam, that is (Ref. 8), 

(j). A {cosh (X. x/L) + cos (X. x/L) - R. [sinh (X. x/L) + sin (X. x/L) ] } 


n 


n 



(i 1 , . • . , n) 

(A. 25) 

< — ^ 
1 

where X. , . . . , X„ 
1 n 

are consecutive roots of the equation 


1 

i 


cos X cosh X - 1 = 0 

(A. 26) 

n 

and 



n 

1 

R. A ( cosh X. - cos X. )/( sinh X. - sin X.) (i = l,..,,n) 

(A. 27) 

. ^ 

i 

then, with p = constant, the integrals in Eqs. (A, 19) through (A. 24) become 


n 

1 

M 

= pL 

(A. 28) 

n 

\ 

d 

= pL^2 = ML/2 

(A, 29) 

] 

n 

I 

= pL^3 - ML^3 

(A. 30) 

i 

1 

1 

= 0 (i = 1, ... , n) 

(A. 31) 
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j 

n 




n 


N. 

1 


0 


(i — 1, •••) 


n) 


(A. 32) 


E.j = 6.^ (i, J = 1, , n) (A. 33) 

where 6.. is the Kronecker delta. 

13 

Equations (A. 15) through (A. 18) thus give way to 


(f«)^ = -UiM + U^UgM + U^d 


(F») ^-UjM-igd-UgUjM 

= - ^3 " 2 ‘>i) - “2 - ^3 Ui d - 2U3 2 % U3^ , 

n 

('^ 3 .j) - -2 1 ") 

i=l 


(A. 34) 
(A. 35) 


(A. 36) 


(A. 37) 


Mounted on B are n^ V-gyro units VI , . . . , Vn^^, a typical one of which , Va , 

is shown in Fig. A-2. The arrangement consists of a gyro casing Go; attached to B , 
and two identical uniform rotors, Da and E a , carried by gimbals whose axes lie in 
A, are parallel to each other, and are "perpendicular" to v in the following sense: 
the mass center of Ga and of each rotor-gimbal assembly are presumed to remain in 
A, with the mass center Pa of Va being a point fixed on v; thus, if b^^ and 
bg^ are unit vectors parallel to A and respectively parallel and perpendicular to v 
at Pa , then the gimbal axes are parallel to b^^ . 
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The gimbals are geared to each other in such a way that the axes of Da and Ea each 
rotate the same amount relative to r , but in opposite senses, as shown. Da and 
Ea are made to rotate with a constant angular speed relative to their respective 
gimbals such that when Yq, == 0 » the angular velocities of Da and Ea in Ga are 
■ ^a -la % -la ’ ^respectively. 

If C denotes the reference frame in which and are fixed, then the angular 
velocity of Ga in C can be expressed as 





(o! 1,..., n^ ) 


(A. 38) 


where Cg ^ x Cg, and Oy/ax)| ^ is the value of by/^x evaluated at x = x , 

^a A e; 

the x-coordinate of Pa. Substitution from Eq. (A. 3) into Eq. (A. 38) and comparison 
with Eq. (A. 7) then gives 


C Ga 
0) 


2 "3+ 1 


(d<^./dx)I Cg 



Hg) 


(A. 39) 


where (d<^./dx) | is the value of d<^ /dx at x = x . Then, noting that the angular 

a 

velocities co of j 3 q, angular velocity 

of C in A are given by 


Ga Da 

CO 


Si 


a 


cos y b, 
'a -la 


U 


3 + n + 


K 

a -2a 


fi sin y c„ 
a 'a -3 


(a 1, ..., n^ ) 
(A. 40) 


81 



°“/“= li, U3 H. „ + „ bj^ - S2„ sinv„ C3 (« = 1, ...,«(,) 


(A. 41) 


AC ,, 

=^303 


(A. 42) 


where 


(q! !)•••) Hq ) (A • 43) 

one can form the angular velocities of Do; and Ea in A by using 

Eqs. (A. 39 ) through (A. 42) and the addition theorem for angular velocities (Ref, 9, 
p. 24), obtaining the expressions 


U 


3 + n + 


A y 

a = 'a 


^CiP^ = - £2 cos y b - U 
— a a -la 3 + n + a -2o! 


n 




i = l 


a 


- S2 sin y 
a a 


C 3 (.a 1 ,,,,, ) 


(A. 44) 


A, ,Eoi _ 


oj = J 2 cos y b. + U„ b 

— a a -1q! 3 + n + q! - 


2 q; 


n 

"3 * 2 "3 4 i 


i = l 


a 


- fi sin y 
a 'a 


C3 (q! = 1, , n^) 

(A. 45) 
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Differentiating Eqs, (A. 44), (A. 45) with respect to t in A then yields the angular 
accelerations and of Do; and Eq; in A, 


A^Do! 


n 


n U„ sinv + U„ U„ + > U„ . (dd)./dx) U„ 

0 ! 3+n+Q; 3 3+n+o; X 3+i ' 3+n+o; 

i _i O' 


i = l 




U„ , , + U„ ^ cos 

3 + n+o; 3 a 


n 


n 

y + / U„ , . (dd>./dx) I n cosy 
'a 3+1 ' ' ' X O' 'a 

01 


i = l 


+ U„ + 7 U„, . (d<A/dx)| - S2 U cosy i^c- 

3 ^ 3+1 ' ^ "x a 3 + n+o; 'o; J -3 

L i=i “ J 


(o; = 1 , . . . , n^) 


-2o; 


(A. 46) 


A^Eo; 


Uo.„. sin 
O' 3+n+o; 


n 

T'«*''3''3+n+a-'-2 “3-h(‘»V*‘)Ix "s+n+a 

.•^1 (X 


i = l 


bi 

-la 


n 


Uo , , + fi cos 

3 + n+Q/ 3 a 


a 


n 

i = l 


i = l 


- J2 U„ . , cos y 

X a 3 + n+o; 'a 
a 


b 


2o; 


Cg (o: = l,...,n^) 


(A. 47) 


Furthermore, Eqs. (A. 44), (A, 45) provide partial angular velocities » 

A. F 

associated with =1, ...,3 + n+n^), these being given by (Ref. 9, 

p. 17) 
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n 


n 


A Da _ A Da _ A Ea _ A Ea _ . 

^1 - ^2 - oji - ^2 - 0 


' — \ 
j 

; 

A Da A Ea 

- ^i?3 “ ®3 




^(a ) 

'n 

■^ 0 )?“. = = (d0/dx)| Cg (j = l,...,n) 

-3+j -J+j J <5 


i 

n 

1 

A Da _ A Ea _ u 

-3+n+a " -3+n+a -2a 


1 

f — 

1 


(A. 48) 

. i 

and, if one assumes that the central inertia ellipsoid of each rotor-gimbal assembly is 
a sphere with associated moment of inertia , then the inertia torques , 

n 

j 

n 

(T*)„ acting respectively on Da and Ea are given by (Ref. i 

"• xUOi 

9, p. 116) 

1 

i 

(or = 1 ng) 

(A. 49) 

n 



'—j 

(T*)j,„ = - V“ (“ = 1 ”0> 

(A. 50) 

1 

! 

A, A. pQf 

The acceleration a and partial velocities v^ ( £ = 1 , 

3 + n + a) of 

' — 1 

Pa in A can be written [see Eqs. (A. 8), (A. 9)] 


n 

I 



, — 1 

i 

i 



^ \ 

! 

j 
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i 


I 



AaPo! 


n n 

i-1 “ i = l “ 


«2 «3 - 4 


)- 


f D. n \ 

"2 " 2 *3 . " «3Ul - «3 2 % ^i'x )SS 

^ 1 = 1 “ 1-1 “/ 


(A. 51) 


= 2i 


(A. 52) 



A Pq; , I , . 1 

-3+j *^j^x -2 


A Po; 
-3+n+a 


0 (or ly***) ) 


(A. 53) 


(A. 54) 


(A. 55) 


(A. 56) 


where (j) | denotes the value of a . at x = x (a; = 1 , . . . , ). 

.1 Oi ’G' 

Finally, if is the mass of Va, and if one assumes that the distance between the 
centers of Da and Ea is small compared with L (a = 1, ... , n^) , then the con- 
tribution (F|)^ from VI,..., Vn^^ to the generalized inertia force associated with 

Ujj (£ = 1, , 3 + n + n^) is given by (Ref. 9, p. 123) 
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h)y = 1 1 ^^ 

a = l '■ 




A Ea 


„ - m “y - 

Ea a -i 


A„Pq! A_Po; 


{i = 1 , . . . , 3 + n + Hp ) (A. 57) 


Substitution from Eqs. (A. 46) through (A. 56) into Eq. (A. 57) then produces 


“g / n n V 

(P*l) =. - 2 “4«1-U3 2 -2^3 2 "3 + l^'x -U2«3-“^a ) 

V a=i \ i=l “ i=l « / 


(A. 58) 


/ n 

(n) = - 2 2 Vi^ii 


Of =1 


i = l 


+ U. X + U. U 
X 3 a 3 

a 


i=i “/ 


(A. 59) 


in). 


*^G r n / n n 

- 2 "“c -2 ‘»j^j'x„ri-“3 2 2 

o; = l L j = l “ \ i = l “ i = l 


a 


-« 2«3 


K* 2 

\ i = l 


+ U„ X + U„ U, 
X 3 o! 3 1 

a 


)1 


n 


G 


Q! =1 


a 


n 


U„ + 7 Uoa.- (d0./dx)| - S2 U„, , cosy 

3 ^ 3 + 1 ' ''x a 3 + n+o! 'a 

i = l 


a 


(A. 60) 
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n. 


G 

V Q! = 1 ' ’ = 1 


n 


i = l 


n \ “g 

i=l ^/ a=l 


J (d<j!)./dx)| 

O' ' ' ' X 


a 


U, 


n 


> Uo , . (d<*./dx) - ^2 U„, , cosy 

/ 3+1 ' ^1 "x a 3 + n+a 'a 

“t 01 

(j = 1 , ... , n) (A. 61) 


i = l 


n 


/ F* . \ = - 2J U„. . + U cosy + X U„ . . (dd./dx) J2 cosy 

\ 3 + n+o!/ 0 ! 3 + n+o: Z a a A 3 + i ' ^i ''x a 

' 'V • 1 a 


a 


i = l 


(a = 1, .... n^) 


(A. 62) 


It is further assumed that the inertia torques acting on Go; (a = 1 , . . . , n^ ) are small 
in comparison with the inertia torques given by Eqs. (A. 49), (A. 50), and, hence, are 
neglected in this analysis. 

Motors attached to Ga exert systems of forces on the gimbals carrying Da and Ea 
equivalent to couples of torques and-T^, respectively, where 

T = T b (a = 1, ..., n^) (A. 63) 

—a a -2a ' ’ ’ G' ' ' 

Thus, the generaiized active force F_ , , associated with U„ . . (a = 1, . . . , n„ ) 

’ “ 3+n+a 3+n+a ' ’ ’ G' 

can be obtained from the expression (Ref. 9, p. 81) 
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F 


3+n+o! 


5 ( 

a = l \ 


A Da rr _ A Ea . „ 

— 3+n+o; —a —3+n+o; —a 


(o; = 1, , np) (A. 64) 


and is found from Eq. (A. 63) and the last of Eq. (A. 48) to be 


=-2T (Q; = l,...,np) 

3 + n+o; a ’ ’ G' 


(A. 65) 


The generalized active forces associated with internal forces in B can be found 
from a strain energy function S and a dissipation function D, given by 


S ^ 


^ 2 2 ^ 
(1/2) J EIO>/9x^) dx 


(A. 66) 


and 




D = (1/2) j cOy/at)^dx 

o 


(A. 67) 


where the flexural rigidity El of B and the structural damping coefficient c of B 
are constants. Substitution from Eq. (A. 3) into Eqs. (A. 66), (A. 67) then gives 


S = (1/2) X) s.- q. q. 

i=l j=l ' J 


(A. 68) 


and 


D = (1/2) i; Z 

i=l j=l ^ J 


(A. 69) 


where 


S.J ^ El 


/ (d^0./dx^) (d^0./dx^) dx (i,j = 1, . . . ,n) (A. 70 


dy = c J 0. 0 dx 

O 


(^>3 l>...,n) 


A. 71) 
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Evaluation of the integrals in Eqs. (A. 70), (A. 71), by reference to Eq. (A. 25), 
produces the relations 


s.. = El xf 6y/(ML^ 


d.. = c L a../M 

1 ] 1 ] 


so that one has, from Eqs. (A. 68), A. 69), and (A. 7), 

n „ „ 


D = [cL/(2M)l U 

i=l 


where 


can be obtained from the expression 


= -9S/3qj - 9D/8U3^j 


which leads to 


-wfq. - cU_. 
1 ^1 3+] 


(i,i = l,...,n) 

(A. 72) 

(i,j = 1, . . . ,n) 

(A. 73) 

.7), 

2 

^i 

(A. 74) 


(A. 75) 

(1 = 1 n) 

(A. 76) 

ted with Ug^. (j = 

1, . . . ,n) 

(j = 1,... ,n) 

(A. 77) 

(j = l,...,n) 

(A. 78) 
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The 3+n+riQ dynamical equations governing the^orce-free) motion of the system 
can now be constructed as follows* (Ref.- 9, p. 177): 



+ 

i% - 

(A. 34) 


(A. 58) 


+ 

1 ! 

(A. 35) 


(A. 59) 

(^3) 

+ 

(^3) = 

\ 3/g 



(A. 36) 


(A. 60) 


+ 


(A. 37) 


(A. 61) 

^^3+n+O! 

+ 

F 

3-t-n+Q; 

(A. 62) 


(A. 65) 


= 0 


(A. 79) 


= 0 


■‘■^3-fj ~ ® (j = l,,..,n) 

(A. 78) 


0 (qj 1 )««.) ^ 


(A. 80) 


(A. 81) 


(A. 82) 


(A. 83) 


A complete description of the motion of the system can be obtained from Eqs. (A. 79) 
through (A. 83) together with the kinematical equations, Eqs. (A. 4) through (A. 7), 
(A. 43) , once initial values of X , Y , 0 , q^, . . . , q^ , 7^, . . . , Tn^ , U . . . , 
have been provided. ^ 


*N umbers in parentheses beneath terms refer to correspondingly numbered 
equations . 
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Appendix B 

LOCAL ROTATION MODES 


In the earlier discussion of Generalized Modes and Force Distributions, reference 
was made, for three-dimensional structures, to "local rotation" mode shapes obtained 
by "stacking" the three-dimensional curls v x of the deformation mode shapes. 

This concept is a direct application of a kinematical theorem due to Helmholtz who 
establishes the following well-known result at the beginning of his paper on vortex motions 
(ref. 10); The most general motion of a sufficiently small element of a deformable body 
can be represented as the composition of 

1. a translation 

2. a rotation 

3. an extension (contraction) in three mutually orthogonal directions 

The proof is based on the multiple Taylor expansion of the relative displacement of two 
neighboring points in terms of their original coordinate differences. In the context of 
linear elastic structures, we will restrict ourselves to differential displacements, ig- 
noring the constant term giving rise to translation. 

We consider 3-d. o.f. finite-element structures in which a point located at r*^ undergoes 
a deformation which displaces it to r^ + u*, as shown in Figure B-1, The deformation 
iT is a time-varying vector field iT = iT(r^, t) , and for finite-element structures, is 
modally decomposed in a set of (usually) orthogonal functions { ^ 

n 
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/ 

0 


Figure B-1. Local structural deformation 

where ( t ) are the modal amplitudes associated with the decomposition in the mode 

shape basis { T } . 

n 

The differential deformation between two points situated at IT and r*^ + dr^ before the 
deformation is then given by 


diT = (viT) • dr" (B.2) 

where V vT is the tensor spatial gradient of the deformation field 1T(7, t ) . This 
tensor may be split into its symmetric and skew-symmetric parts (using Gibbs' dyadic 
notation); 


VtT = 


S : Strain Tensor 
(Symmetric) 


+ 


I (v7-lTv) 

' ^ ^ 

I ( viT xiT) X 1 

(Skew-Symmetric) 


(B.3) 


where 1 is the identity dyadic. If we define; 
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60 = I (V X u) 


(B.4) 


, the differential deformation Eq. (B. 2) can be written as 

diT = l)* d?* 

= S • diT + ^ X dr^ 


(B.5) 


and 66 = 66 (r, t ) is thus the instantaneous axis of rotation of a local material 
neighborhood of a point of the structure. 

To obtain the "local rotation" modes, we substitute Eq. (B. 1) into Eq. (B.4) to obtain 


where 


60 


= = 2 V‘> 


n 


= 2 


qn<t) 50^(r) 


n 


a\(7) . i (v X 


(B.6) 


(B.7) 


is the nth modal rotation coefficient , i.e. , the contribution of the nth mode to the local 
instantaneous axis of rotation at a point iT. The corresponding nth column of the local 
rotation modal matrix is then obtained by "stacking" columnwise the three components 
of the Gibbs vector 60^ according to the nodal sequence of structural stations employed 
in the finite -element decomposition initially used for the • 
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MODAL TILT COEFFICIENTS FOR SEGMENTS DEFINED BY n POINTS 


A special case of the modal rotation coefficients occurs for flat, plate-like structures 
for which the instantaneous local rotation axis remains essentially tangent to the struc- 
ture. This arises, for instance, when high-frequency shearing modes are ignored in 
a given structural model. This situation, in itself, does not justify a special formula 
simply because ^ = ^(^x^) now only has two nonzero components in the local 

frame. If, however, for a structure of this kind, several relatively well separated 
structural nodal points are attached to a rigid plate (supporting, for example, a set of 
gyrodampers) , then the "small neighborhood" assumption required for Helmholtz’s 
theorem is no longer valid. 

To associate a local rotation (or tilt) axis with a finite rigid plate interconnecting a 
number of n distinct structural nodal points, one can, of course, resort to 6-d. o. f. 
finite-element models. The difficulty with this procedure comes from the numerical 
methodology employed in synthesizing (or "compressing") model data into the multibody 
nonlinear dynamics codes, which incorporate flexibility with the usual 3-d. o.f. finite- 
element models. A very good approximation can be obtained , for 3-d. o.f. finite 
element models, by assuming that the local deformation will move the rigid plate into 
an orientation which is assumed to coincide with the least squares best-fit plane for 
the n nodal points, deforming individually without the constraint of being bolted to a 
rigid plate. 

The algorithm for this procedure is displayed in Figure B-2, and numerical compari- 
sons with 6-d. o.f. finite -element models have established the validity of the procedure 
for the special "flat" structures considered. 
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3 xn MATRIX: X = f?, , . . . ,7 1 (n POINTS OF DEFORMED STRUCTURE) 

I n. 


f3 xn MATRIX 


IFORJTHMODE; I 


= p' (JTH modal matrix at those n POINTS) 


\ ^CENTROI^/^% 




LEAST SQUARES 
BEST FIT PLANE _ 
WITH NORMAL 7 


DEFORMED 

STRUCTURE 

SEGMENT 


X^ s (XX^) ' X (PSEUDOINVERSE OF X , 3 xn) 

7 = xM'-i-.-wir 

q.(0: AMPLITUDE OF JTH MODE 
• TILT DEFINED BY: 


= -— p- X lx"' IdXl^T] 

I7|2 171 2 I ' 

WHERE DEFORMATIONS (dXj = X!.q.(t)(*.j 


DEFINE "MODAL TILT COEFFICIENTS" de. BY 

, 

6J, = - _L_ X [x^ ^-*>.)^7] FOR JTH MODE 


SO THAT TOTAL TILT IS GIVEN BY: 


69(0 = 2 RjW 


Figure B-2. Modal tilt coefficients for segments defined by n points . 


95 





Appendix C 

GYRODAMPER SCALING EXAMPLE 



The following example illustrates the methodology for scaling a structural control 
system using gyrodampers. It is assumed that the structure is given, and that 
an existing CMG device, with acceptable performance, needs to be scaled to 
match the structure. In this instance, the structure is chosen to be a free- 
free aluminum beam with the following characteristics: 

M = 154.6 kg 
L = 6.9 m 
I = 4.375 10"®m‘^ 

E = 7.2 10^® Pa 

P = 2760 kg/m^ 

(K„= (EI)-1/2 = 1.782 10-3) 

s 


The V-gyrodamper(s) will be located at one (both) ends of the beam. The scaling 
procedure is as follows: 

Step 1 . The values of the rotational mode shapes at the beam-ends are computed 

for the first few modes of interest (either analytically or by finite-element models) . 

- 1/2 -1 

In this case, the first rotational mode shape is 0.1098 rad kg m with 

a frequency fj^ = 8.8 Hz. 

Step 2 . The damping gain D is calculated from the desired damping ratio J 
using formula (16). In this case, a 10% damping ratio is desired in the first 
mode, leading to a gain 

2^0) 

D = ^ = 916.55 Nms 
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A full complex eigenanalysis of the system provides then then damping of the 
various modes. If not satisfactory, D is modified accordingly. 

Step 3 . An existing gyro design is selected. In this case, the CMG chosen is 
the Bendix Model No. MAS- 100-1 shown in Table C-1 below. Its characteristics 
are: 


m = 
h 
n 
J 


*^max 



17.3 kg 
6.78 Nms 
8000 rpm 
0.01 kg m^ 
20 rad /sec 

16.3 Nm 
320 watts 


From Eq. (54), the following GSP are computed for this class of designs: 


K 

K 

K 

K 

K 


H 

R 

J 

G 

P 


0.059 

214 

37,450 

1.523 

29.91 


^ GSP Values 


Step 4 . The mass ratio fx is chosen on the basis of an acceptable weight for the 
control system, e.g., p = 0.1. Since M = from Eq. (64), the number N of 

V-CMG pairs can be increased provided smaller gyros are used such that the 
product Nm remains constant. The best value for N is determined then by 
examining the four limits given by Eq. (67). The corresponding composite graph 
of these inequalities is shown in Figure C-1, where the maximum vibration 
amplitude 0^^ is plotted versus the number N on a log-log scale. As can be 
seen in this particular example, the optimum choice is N = 2 V-CMG pairs. 
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Table C-1 

BENDIX CONTROL MOMENT GYROS 


MODEL NO: 

MA-1000 

MA-2300 

MA-500 AC 

MA-500 DC 

MA-2000 

MA-5-100-1 

TYPE 

DOUBLE GIMBAL 

DOUBLE GIMBAL 

SINGLE GIMBAL 

SINGLE GIMBAL 

DOUBLE GIMBAL 

SINGLE GIMBAL 

PROGRAM 

CMG R and O 

Skylab 

CMG Test 

CMG R and 0 

ADV.CMG 

Pointer Control 

AGENCY 

NASA LRC 

NASA MSFC 

Air Force 

In-House 

NASA MSFC 

LMSC 

UNITS BUILT 

4 

12 

3 

5 

1 

5 

STATUS 

Evaluation at 
NASA 

Qualified 

Partially 

Qualified 

Partially 

Qualified 

Evaluation at 
NASA 

Evaluation at 
LMSC 

MOMENTUM 

Utlb-»ec) 

tooo 

2300 

250* 750 

250 - 750 

1000-3000 

5 

OUTPUT TORQUE 
Max ft'ib 

175 

122 

500 

500 

. 175 

100 

PIVOT TORQUE 
Max ft'lb 

175 

122 

90 

90 

175 

12 

DECREES OF 
FREEDOM 

2 

2 

1 

1 

2 

1 

GIM0AL FREEDOM 
(deg) 

Unlimited — 
Slip Ring» 

±80'*±175*' 

+ 170° 

+ 170° 

Unlimited — 
Slip Rings 

Unlimited — 
Slip Rings 

GIMBAL RATE 
Max (dcg/scc) 

10 

4.7 

57.3 

57.3 

5,30 

1 146 

WEIGHT 

Obs) 

230 

418 

145 

155 

558 

38 

APPROX. 

ENVELOPE 

39" Oiam. 
Sphere 

4 1 " Oiam. 
Sphere 

(Cylinder 20" 

Oiam. X 32" Long) 

4‘4" Oiam. 
Sphere 

Cylinder 10" x 
10" Long 

ROTOR SPEED 
(RPM) 

11,400 

9,000 

7,850 

7.850 

4.000 to 

12.000 

8.000 

SPIN UP 
TIME (hr$) 

<^4 hours 

14 hours 

<C6 hours 

<C 2.5 hours 

2.5 hours 
(to 8.000 RPM) 

<0.1 


J _..J __ 


J 



J 




J 


- J 


J 


j 



] 

1 






step 5 . The characteristics of the required V- gyros are then obtainod from 
their individual mass, i.e., 


= 7.73 kg 

which is used in Eq. (54) with the GSP values given above. The geometric 
reduction factor is then ^ 7. 73/17. 3 = 0.764. 

Step 6 . If the desired is not high enough, the ratio /u can be increased 
(which will push all the limits upward). Otherwise, by examining the graph, 
manufacturer design modifications can be sought in order to meet desired per- 
formance. Eq. (67) could be used at that point to determine parametrically 
which gyro characteristics need to be improved. 
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